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Abstract 

We calculate the vrSA and vrES coupling constants in light cone QCD sum rules 
for the structure o'apJbP^'Q^ ■ A comparison of our results on these coupling constants 
with prediction of the SU(3) symmetry is presented. 
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1 Introduction 



Many strong interaction-processes involve meson-baryon coupling constants as the main 
ingredient. The determination of these fundamental quantities requires information about 
the physics at large distance. In other words, for a reliable determination of these parame- 
ters we need some nonperturbative approach. Among all nonperturbative approaches, QCD 
sum rules is one of the most powerful method in studying the properties of hadrons. 
This method is based on the short distance OPE of vacuum-vacuum correlation function 
in terms of condensates. For the processes involving light mesons vr, K or p, there is an 
alternative method to the traditional QCD sum rules, namely, light cone QCD sum rules 
0. In this approach the expansion of the vacuum-meson correlator is performed near the 
light cone in terms of the meson wave functions. The meson wave functions are defined by 
the matrix elements of non-local composite operators sandwiched between the meson and 
vacuum states and classified by their twists, rather than dimensions of the operators, as is 
the case in the traditional sum rules. Many applications of light-cone QCD sum rules can 
be found in [Q-Jlll and references therein. 



In this work we use light cone QCD sum rules approach for determination of the coupling 
constants of the pion to the lowest states of the baryon octet S and A, q-^y.^ and S'ttSE- 
Note that these coupling constants were investigated in framework of the QCD sum rules 
based on pion-to-vacuum matrix element in the leading order of the pion momentum q 



for the structure ^75 in where q is the pion momentum. The results of this work 



are currently under debate in literature (see discussions in and 0]). Moreover in |jT5 
and ||T6| it was pointed out that there is coupling scheme dependence for the structures 75, 
^75, i.e., dependence on the pseudoscalar or pseudovector forms of the effective interaction 
Lagrangian of pion with hadrons have been used, while the structure (T^i/75 is shown to 
be independent of any coupling schemes. For this reason, in present work we choose the 
structure cr^y^^p^q^ , where p and q are the A (E) and the pion momenta, respectively. 
It should be noted that the sum rules for the CF^y'-^^p^q" structure was derived in [1^ in 
investigation of (^p^;^ coupling constant. 

The paper is organized as follows. In Section 2 we derive sum rules for the pion-baryon 
coupling constants (7,rEA and q-kt^t. for the structure (J^ulbP'^Q'^ ■ Section 3 is devoted to the 
numerical analysis of the sum rules for gnT.A and s^ttss and discussion. 

2 Formulation of the pion— baryon sum rule for the 

^^SA and ^^ss 

According to the main philosophy of the QCD sum rules, a quantitative estimation for gnj:A 
and gn-Es couplings can be obtained by matching the representations of a suitable correlator 
in terms of hadronic (physical part) and quark-gluon language (theoretical part). For this 
purpose we consider the following two-point correlator function with pion 

U{p,q) = f rf^xe^^^(7r(g)|T[r/y(x)r/s+(0)]|0) , (1) 



where p and rjY are the four-momentum of the hyperon (in our case A° or S°) and its 
interpolation current, respectively, ?7s+ is the interpolating current of and q is the pion 
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four-momentum. The interpolating currents for A°, and S+ are 



r]YO = a eabc 



VT,+ = ^abc (u^C-f^Ub) l5l''Sc , (2) 
where s, u and d are strange, up and down quark fields, the upper(lower) sign corresponds 



to A°(S°) and a = y^2/3 for hP and ^2 for S°, respectively, a, 6, c are the color indices, C 
is the charge conjugation operator. Saturating correlator (1) with the Y (= A° or TP) and 
states in the phenomeno logical part, we get 

n = ('(.)yf> (s> + ,)|,fe.|o>(o|,>|r(ri) ^ ^^^^^^ 

The matrix elements in Eq. (3) are defined in the following way 

{Q\r^Y{x)\Y{p)) = Xyu{p) , 

'7r(g)r(p)|S+(p + g)) = -gYi^+n- HphMP + ■ (4) 



Substituting Eq. (4) in Eq. (3), and choosing the structure i (Jai3p"q^'y5 for the physical 
part of Eq. (1), we get 

Yjphys ^ ——- rf/'^^'^fo^ + higher resonances . (5) 

[p^ — [{p + qY — 

Let us now consider the theoretical part of the correlator (1). From this correlator we 
have (we present only the terms which give contributions to the above-mentioned Lorentz 
structure) 

n = -ay'd^a;e'P"{-757M757^7pC5^C-S^5^7p75(7r|n757^rf|0) 
±[ - 757/.'57pC(757v)^C~Sm'5'7p75 (tt \u-f5l^d\ 0) 

+^l5lf^S-fp CalpC'^-ipS^-ip-i^ (tt \uaapd\ 0) ] } , (6) 



where upper (lower) sign corresponds to A(S) case and a = y 2/3 for A while a = \pl for 
E. Here S and iS* are the propagators containing both perturbative and nonperturbative 
contributions, respectively. Here we present the explicit form of iS^{x) 

■osl \ ^ i "^'^ 1 ^ I- \ (^ ^'^s , 
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mo(ss) 1 
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where is the mass of the strange quark and Gap is the gluon field strength tensor. 
The form of S can be obtained from Eq. (7) by making the replacements (ss) — > {qq) 
and m<( ^ 0. From Eq. (6) we observe that, in calculation of the correlator function in 
QCD, the matrix element of the nonlocal operators between the vacuum and pion states 
are needed. These matrix elements define the two particle pion wave functions and up to 
twist four they can be written as 







^/tt^m / due''"^''[ip^{u)+x^gi{u)] 



7r(g) (i(x)cr„/375M(0) 







qx 



due''"^''g2{u) , 



6(m„ + rrid) Jo 



due''"^''ip^{u) . (8) 



Using Eqs. (6), (7) and (8) we get the following result for theoretical part (for the structure 

Yjtheor 



+ 
+ 
where 



+ (-X±a) 



+ 



TTln 



(-A±a) + (1±1) 



230471^x2 



1447r2x2 727r2a;2 



du e 



luqx 



27r%4 



c/Me*"'?"[(7i(M) + G2(M)] 



(9) 



A 
a 



92{v)dv 

{qq) - (ss) , 
{qq) + (ss) . 



Our next task is to perform integration over x and perform double Borel transformation 
in Eq. (9) with respect to the variables p"^ and {p + g)^, in order to get an answer for the 
theoretical part of the sum rules. As an example, let us demonstrate on one of the terms in 
Eq. (9) how integration over x and double Borel transformation can be carried. Consider 
the following term 



du(pT,{u) / d'^x 



X^ 



-laapXaqp . 



In performing the x integration, we will make use of the formula 

d^x 



(X2)" 



z(-l)"2 
V{n - l)T{n) 



(pT-'M-r)+^n-2 , {n>2) 
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where Vn-2 is a polynomial of power n — 2. However this polynomial is inessential, since it 
vanishes after double Borel transformation. Hence, disregarding this polynomial we have 



d X- 



oi(p+qu)x 



f Pa 

= Sy du(p^{u)^iaapqi3 , 

where P = p + qu. The last step in this calculation is performing double Borel transforma- 
tion over the variables and {p + q^ to the expression 

/ du^^{u)- — -2 . 

(p + qu) 

Rewriting {p + qu)"^ = p^u + u{j) + g)^ (here the pion mass is neglected) and using the 
exponential representation for the denominator, 

Jo 



we have 



j du(pT^{u) 



\p^u + u{p + qj" 
1 



(p + quY 



The double Borel transformation over the variable p^ and {p + q)^ is done with the help of 
the following general formula 



Tin) 



where 



and Uq 



{M^f--8{u-uo) , 
Ml 



Ml + M| ^ Ml + M| ' 

and Ml and M| are the Borel parameters. Using this expression and performing the 
integrations over the variables a and m, we finally get 

du(pT,{u) I d'^x- 



x^ 



All other terms can be calculated similarly and for the theoretical part it follows from Eq. 

(9) that 

Y^theor 



+ [gi{uo) + G2{uo)] 



(1 ± l)£jMVi(«o/M^) - ^(-A ± a)M'fo{so/M') + ^(-A ± a) 



(1 ^ 1)^MVi(WM2) + (1 ± + (1 ^ i)LmJss) ± -m,(7 

^ ^ ^247r2 ■'^^ ^ ^ ^2304 v ^ y-^g s\ / g ^ 



;i ± l)\msM^fo{so/M^) + ^ (-A ± a) 



(10) 



where the function 

2 ^ {sq/M^Y 



is the factor used to subtract the continuum, which is modeled by the dispersion integral 
in the region Si, S2 > Sq, Sq being the continuum threshold (obviously the continuum 
thresholds for the A and S channels are different). Since masses of A and S are very close 
to each other, we can choose Mf and M| to be equal to each other, i.e., Mf = M| = 2M^, 
from which it follows that Uq = 1/2. 

Performing double Borel transformation over the variables p"^ and (p+g)^ in the physical 
part (5) and then equating the the obtained result to Eq. (10), we get the sum rules for 
S'ttAe and gnST: coupling constants 

^yE+.-AyAs+=e-'/^'^*''-^ (11) 

where m ~ rrix ^ rriY.- 

From Eq. (11) it follows that in determining the strong coupling constants QnAT. and 
S'ttSS the experimentally undetermined residues Aa and Ae need to be eliminated from sum 
rules. The residues Aa and As are determined from corresponding mass sum rules for the 
A and S hyperons |]18|, as follows 



+6MVo(4/^') + ^«'(3 + 47) , (12) 
\X^fe-"^y'''327r^ = /M') - 2am,{l + ^)M^ fo{s^ /M^) 

+bM^fo{s^/M^) + ^a^ , (13) 



where 



a = -27r (qq) , 

b = ~ 0.012 Gel^\ 

and the functions fo{x), fi{x) are presented just after Eq. (10). The ratio of the Eqs. (11) 
(12) and (13) gives 



Yltheor 



yttheor 

S'ttSS — e —2 — . [io) 



The main reason why we consider the above-mentioned ratio rather than the individual 
sum rules themselves (i.e., first determine Aa and As independently from Eqs. (12) and 
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(13) and substitute their obtained values in Eq. (11)) is that the sum rules obtained from 
these ratios are more stable as is similar to the baryon mass sum rules case. In addition 
to that the uncertainties coming from various parameters such as quark condensate, m^, 
continuum threshold sq and Borel parameter, are reduced. 



3 Numerical analysis 



Now we are ready to perform the numerical analysis. The main nonperturbative input 
parameters in the sum rules (11) are the pion wave functions. In our calculations we have 
used the set of wave functions proposed in . The explicit expressions of the wave functions 
are 



1 + a2{^i)Cl''^{2u - 1) + a4(/i)Cf 2(2n - 1 

15 



hiu-u? - 1 



Quu 



+ 21(m - u f - U{u -u^ + l 

8 L 

uu{2 + 13uu) 



G2(M,/i) 



+ 10m^ lnu(2-3u+ ^u^^ + lOu^ In m (^2 - 3m + ^u^ 

5 _2 / \ — 2 2 

-6 {fiju u , 



(16) 



where u = 1 — u, 

CT and CT 

are the Gegenbauer polynomials defined as 



3 r 
2 

1) = H 



5(2m- 1)2 + 1 



21(2m-1)^-14(2m- 1)2 + 1 



;i7) 



and a2(/i = 0.5 GeV^ = 2/3, a^^fi = 0.5 GeV) = 0.43. The parameters 6'^{fi 
0.2 GeV^ |211 and eifi = 1 GeV) = 0.5 §. Furthermore U = 0-132 GeV, /i^(/i 



1 GeV) 
1 GeV) 



1.65, {ss) = 0.8{qq) and {qq)\^,=l ceV = -(0.243)^ GeV^, sq = 4 - = (3.0 ±0.2) GeV^. 
Moreover all further calculations are performed at m = mq = 1/2. 

Having fixed the input parameters, one must find the range of values of over which 
the sum rule is reliable. The lowest possible value of is determined by the requirement 
that the terms proportional to the highest inverse power of the Borel parameters stay 
reasonable small. The upper bound of is determined by demanding that the continuum 
contribution is not too large. The interval of which satisfies both conditions is 1 GeV^ < 
< 2.5 GeV^. The analysis of the sum rules (14) and (15) shows that the best stability 
is achieved in the region of M^, 1.4 GeV'^ < < 1.8 GeV'^. This leads to the following 
result for the coupling constants 



5'7rAS 



5.3 ± 1.8 
12.5 ±4.5 
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in which the errors coming from the quark condensate (which varies in the region —(0.24 GeVY 
and —(0.26 GeVY), rn^ parameter (which varies in the region 0.6 GeV^ and 1.4 GeV"^), 
variation of the Borel parameter and change of the continuum threshold sq are all taken 
into consideration. Our calculation shows that the main error comes from uncertainties 
of the quark condensate. The central values of the coupling constants are obtained at 
ml = 0.8 GeV^, (qq) = -(0.243 GeVf and Sq = 3.0 GeV\ 

At this point we would like to stress that the above-mentioned strong coupling constant 
have been analyzed using the individual sum rules themselves by first determining Aa and 
As independently from Eqs. (12) and (13) and substituting their obtained values in Eq. 
(11). The results predicted by this approach are close to the ones presented in Eq. (18), 
however our calculations show that the ratio sum rules prediction is more stable and reliable. 

Here it should be noted that since the phase of the coupling constants can not be 
predicted by the sum rules, we take into consideration the magnitudes of these coupling 
constants to be able to compare them with the predictions of other approaches. 

Let us now compare our results on the coupling constants gnAE and gnj^T: with SU(3) 
symmetry prediction. As is known, SU(3) symmetry predicts 

2 

CttAS = y^(l — a) GjVATTr 5 

GttEE = 2aGjVAr7r , (19) 

where a = F/{F + D) (see for example p2|)- Exact SU(3) symmetric analysis of pion- 
baryon coupling gives F/D ~ 0.58 [^] (exact SU(6) symmetry predicts F/D = 2/3). It 
follows from Eqs. (19) that 

fl = ^ = ^«l. (20) 

while our analysis yields i? ~ 2. It follows from a comparison of these results that SU(3) 
symmetry is broken significantly. The pion-baryon couplings predicted in fl^ are not 
presented here since we have already noted that the results of this work is currently under 
debate in literature. 

In conclusion we have calculated the strong coupling constants of pion with A and S 
hyperons and found out that our results differ significantly from that of the SU(3) symmetry 
prediction. 
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